The physics of jets in FR I radio galaxies 
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Abstract 

We model jets in low-luminosity (FRI) radio galaxies as intrinsically symmetrical, 
axisymmetric, decelerating relativistic flows with transverse velocity gradients. This 
allows us to derive velocity fields and the three-dimensional distributions of magnetic- 
field ordering and rest-frame emissivity. A conservation-law analysis, combining the 
kinematic model with X-ray observations of the surrounding IGM, gives the profiles 
of internal density, pressure, Mach number and entrainment rate along the jets. We 
summarize our recently-published results on 3C 31 and outline new work on other 
sources and adiabatic jet models. 

Key words: galaxies: jets, radio continuum:galaxies. X-rays: galaxies, magnetic 
fields, polarization 
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1 Introduction 



This paper is a progress report on a 
project whose aim is to derive quan- 
titative estimates for the physical 
parameters of jets. We have devel- 
oped techniques to derive the three- 
dimensional distributions of velocity, 
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rest-frame emissivity and magnetic- 
field structure and hence to deduce 
the jet dynamics via a conservation- 
law approach. Our fundamental as- 
sumption is that jets may be approx- 
imated as intrinsically symmetrical, 
time-stationary, axisymmetric rel- 
ativistic flows. We model their ob- 
served radio synchrotron emission in 
total intensity and linear polariza- 
tion, using the observed differences 
between approaching and receding 
jets to constrain velocity, emissiv- 
ity and field. We then combine this 
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model of jet kinematics with a de- 
scription of the surrounding IGM 
and use conservation of energy, mo- 
mentum and particles to estimate 
the internal pressure and density. 
For this technique to work, we need 
sources with two-sided but asymmet- 
rical and straight radio jets. We then 
make deep observations with many 
resolution elements to derive total 
intensity and linear polarization, cor- 
rected for any Faraday rotation. We 
also require X-ray observations of the 
hot gas surrounding the radio source 
from which the external pressure 
and temperature may be derived. At 
present, these requirements are met 
only by VLA and Chandra observa- 
tions of nearby FRI radio galaxies. 



2 3C31 



The first source to which we have ap- 
plied these techniques in full is 3C 31: 
the kinematic model. X-ray observa- 
tions and conservation-law analysis 
are described by p^aing fc Bridle ] 
(|2002a|) , Pardcastle et al. | (g002D 
and |Laing fc Bridle | ( p002b| ), respec- 
tively. The kinematic model requires 
the jets to heatOr^i 52° to the line of 
sight. The jets may be divided into 
three distinct sections by geometry, 
velocity and emissivity variation: a 
well-collimated inner region, a flar- 
ing region in which the jets widen 
rapidly and then recoUimate and a 
conical outer region. Both decelera- 
tion (from /? = f/c ~ 0.8 where the 
jet flares to /3 ~ 0.25 at the end of 
the modelled region) and transverse 
velocity gradients are inferred. The 



toroidal magnetic field component is 
larger than the longitudinal compo- 
nent almost everywhere and the ra- 
dial component is small except at the 
edge of the jet in the flaring region, 
where the field is roughly isotropic, 
perhaps clS 8b result of turbulent en- 
trainment. The beginning of the flar- 
ing region marks a discontinuity in 
emissivity and, probably, velocity. 

The conservation-law analysis shows 
that the jets must be extremely light 
(p 10~^^kgm~^). The entrainment 
rate has a local maximum where the 
expansion is fastest; thereafter, it in- 
creases smoothly and monotonically. 
The required entrainment could be 
provided by stellar mass loss close 
to the nucleus, but interaction with 
the external medium is required at 
larger distances. The jet is overpres- 
sured with respect to the surround- 
ing medium where it flares, but the 
outer region is likely to be in pres- 
sure equilibrium. The composition 
of the jet is not determined uniquely 
by this analysis, but an electron- 
positron jet with entrained thermal 
matter would be consistent with all 
of the available evidence. 



3 New results 



We have applied simplified kine- 
matic models to describe less de- 
tailed observations of a complete 
sample of FRI radio galaxies from 
the B2 sample ([Laing et al. 



Two of these, 0326+39 and 1553+24, 
have been reobserved with the VLA 
at 8.4 GHz and modelled in detail 
(Canvin & Laing, in preparation). 



2 



(a) 0326+39 




lOarcsec p = 1 



Fig. 1. Contours of total intensity superimposed on vectors with lengths proportional 
to the degree of polarization, p, and directions representing the apparent magnetic 
field, (a) VLA observations; (b) model. The resolution is O.Sarcsec FWHM and 
the contour levels are 1, 2, 4, 8, 16, ... x 25 /iJy / beam area. The angular and 
polarization vector scales are shown by the labelled bars. 



Good fits were again obtained: we 
show a comparison between model 
and data for 0326+39 in Fig. 1. The 
basic picture of jet deceleration with 
tranverse velocity gradients holds for 



both sources, but there are interest- 



ing differences in the field structures: 
1553+24 (like 3C31) has a dominant 
toroidal component but the outer 
part of 0326+29 (Fig. 1) has roughly 
equal radial and toroidal compo- 
nents, but no longitudinal field (cf. 
Laing | , |T980| ). 



Thus far, we have assumed that 
emissivity and field structure are in- 
dependent of the velocity field. This 
is appropriate if we aim to deduce 
the internal parameters without im- 
posing preconceptions about the un- 



derlying physics, but to make further 
progress, we need to make additional 
assumptions. The simplest approx- 
imation is that the jets are "adia- 
batic" in the sense defined by [Burch 
(|1979|) : the relativistic particles lose 
energy only by the adiabatic process 
and the magnetic field is convected 
passively with the flow. [Baum et al 



(|1997|) derived analytical relations 
for the surface-brightness of a rela- 
tivistic jet with no transverse velocity 
gradient and either a transverse or 
longitudinal magnetic fleld. We have 
generalized their approach to include 
shear and arbitrary initial field ge- 
ometry, using a formalism based on 
that of liVlatthews & Scheuer 



In Fig. 2, we show the results of fit- 
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Fig. 2. Vectors with magnitudes pro- 
portional to the degree of polarization 
and directions of the apparent magnetic 
field superimposed on contours of total 
intensity at a resolution of 0.75 arcsec 
for the outer regions of the jets in 3C 31. 
The plots cover the range 10 - 27 arcsec 
on either side of the nucleus. Right: 
main jet; left: counter-jet. From the 
top: observations; free Gaussian model 
from Laing &; Bridle ( 2002a| ); adiabatic 
model with initial conditions set at the 
outer boundary and optimized velocity 
and angle to the line of sight; adiabatic 
model with geometry and velocity from 
free model. 



ting adiabatic models to the outer 
regions of the jets in 3C31, with 
initial conditions set as profiles of 
emissivity and field-ordering across 
the jet at the start of the region. 
The models fit reasonably well, but 
cannot accurately describe the ob- 
served polarization at the jet edges. 
A likely possibility is that the veloc- 
ity field is more complicated than 
the simple laminar flow we assume, 
and that some turbulent component 
is present, leading to changes in field 
ordering and strength which are not 
described by our model. In contrast, 
adiabatic models fail completely in 
the inner and flaring regions. This 
should not come as a surprise: the 
X-ray emission detected in these re- 
gions by pardcastle et al. 1 (|2002|) is 
most likely to be synchrotron radi- 
ation, requiring significant particle 
acceleration. 
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